Scaling Behavior in Turbulence is Doubly Anomalous 
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It is shown that the description of anomalous scaling in 
turbulent systems requires the simultaneous use of two nor- 
malization scales. This phenomenon stems from the existence 
of two independent (infinite) sets of anomalous scaling expo- 
nents that appear in leading order, one set due to infrared 
anomalies, and the other due to ultraviolet anomalies. To ex- 
pose this clearly we introduce here a set of local fields whose 
correlation functions depend simultaneously on the the two 
sets of exponents. Thus the Kolmogorov picture of "inertial 
range" scaling is shown to fail because of anomalies that are 
sensitive to the two ends of this range. 
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Anomalous multi-scaling in turbulence is usually dis- 
cussed [Q,^ in terms of the simultaneous structure func- 
tions of velocity differences across a scale R: 



SniR) = (|u(r + R) - u(r)r) ~ (e-i?)"/^ , - 



<5„ 



(1) 



where (...) stands for a suitably defined ensemble av- 
erage, e is the mean energy flux per unit time per unit 
mass, and Sn is the deviation of the scaling exponent C.n of 
the structure function from the Kolmogorov 1941 (K41) 
prediction Cn = n/S — Sn- Since K41 follows from di- 
mensional analysis deviations require a renormaliza- 
tion scale, and it is accepted |^,^,^-|| that in 5„(i?) it is 
the outer scale of turbulence L that serves this purpose. 
The same renormalization scale appears in the correla- 
tion function of the energy dissipation rate e(r, t) ( which 
is roughly i/|Vu(r, with ly the kinematic viscosity) [Q: 



KUR) = Mr + R)- e-)(e(r) - e)) ^ - 



(2) 



where /i is known as the "intermittency exponent" [Q. 
The appearance of the outer renormalization scale in 
these quantities has been correctly interpreted as a fail- 
ure of the K41 basic assumption of inertial range scaling. 
The aim of this Letter is to discuss infinite sets of lo- 
cal turbulent fields whose correlation functions require 
two simultaneous renormalization scales, L and rj where 
r] is the inner (viscous) scale. These correlation func- 
tions demonstrate that K41 fails doubly, once because 
of infrared and once due to ultraviolet anomalies. This 
double anomaly results, in addition to an infinite set of 
multiscaling exponents Cm with a second infinite set of 
exponents that are denoted here as Pi . The phenomenon 
occurs in a similar fashion in scalar turbulent advection 



and in Navier-Stokes turbulence. Since the development 
of the ideas is simpler in the case of scalar fields we will 
present them in the context of scalar advection and gen- 
eralize later to turbulent vector fields. In all cases the 
considerations are based on the fusion |^ of two or more 
near-by coordinates as is shown next. 

First we generate local fields that originate from the 
fusion of two points. Consider for that a turbulent scalar 
field r(r, t) and the product of two such fields at two 
adjacent points 



vI'(p,r) = T(r + p/2)T(r-p/2) . 



(3) 



It is advantageous to represent this field as a multipole 
expansion 



I 



*,(p,r)= ^ Yi„,ip) J ^iptr)YU^)d^ 

m— — l 



(4) 



(5) 



Here p = pj p and ^ are unit vectors. The orthonormal 
spherical harmonics Yim{p) are the eigenfunctions of the 
angular momentum operator L = — ipx V which depends 
only on the direction of p.: 



PYi„,{p)^l{l + l)YUp) 



(6) 



Next we wish to represent 5';(p, r) in terms of (infinitely 
many) local fields which depend on r only. To this aim 
we expand ^'(p, r) in a Taylor series in p. This turns 
Eq.(|) to 



7n— — l 
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n=0 



V) 



'\2r, 



(2n)! 



-*(p',r) 



p'=0 



(7) 



Here and below the operator V'^ — d/dp'^. Note that we 
have only even n orders since our field \E'(p, r) is even in 
p. Performing the angular integrations we end up with 



*,(p,r) J2 



n=l/2 



(2n)! 



(p2v'2)"-'/2$(p',r) 



p'=0 



(8) 



The operators Ai in this equations are related to the ir- 
reducible representation of the SO (3) group in a manner 
that will clarify soon. For the first values of I they are 
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A, = {p- vr - ^p'^'Hp ■ + |p'v'4 . (9) 

The coefficients are determined by requiring orthogonal- 
ity of Ai to aU All with Z' < Z in the sense of integrating 
AiAii over the direction of p. In fact the coefficients of 
the (p • V')P term in Ai are those of x'' in the Legendre 
polynomial Pi{x) On the other hand, the coefficients 
Q2n,i in (jsj) are determined by requiring that the 
will agree with the original Taylor expansion (^) which 
has only orders of (p ■ V). This means for example that 
o-p.p = 1 for any p, 02,0 — 1/3, a4_o = —3/35, etc. Equa- 
tion M) can be rearranged in the form 



EO-l+2p,2p 2p T 



r aia2...ai / 
^l,p \ 



p=0 



(10) 



Here we have introduced the tensorial local fields Li^p 
and the local differential operators Di(r) which for the 
first values of I are 

Do = l 



^4 
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{dapS^S + <5a7<5/35 + SaS^Pa-f) ■ (H) 



Here V^j = d/dpa- The readers familiar with the repre- 
sentations of Lie groups recognize immediately that our 
local fields L/_o are nothing but the 21 + 1 rank irreducible 
representations of the SO (3) group The tensor fields 
thus obtained are symmetric to any pairwise exchange of 
indices. We will propose now that these gradient fields 
have ry-related anomalous scaling which is governed by 
a set of anomalous exponents /?;. Autocorrelation func- 
tions of these fields, and correlation functions of these 
fields together with field differences across a scale R de- 
pend also on R/L with exponents determined by the set 

Denote the correlation function of the tensorial field 
L/_p with 2n — 2 scalar T-fields as C2„.;(r, ra, . . . r2„): 

C2n,i,p(r,r3,...r2„) = (L/,p(r)T(r3)...r(r2„)) . (12) 

Note that in this correlation function p does not appear. 
However it is related to the standard 2?i-point correlation 
function in which two coordinates (say ri and r2) are 
separated by a small distance p. By definition 



p p 

■^2n(r+-,r--,r3, 



■r2n) 



To connect the functions (|T^) and (|T^) we represent 

as a multipole decomposition = Si^o-^2n,i- Using 

(S) we have 



n- / P P 

^2«,i(r + -,r - -,r3 . . . ra,, 
{^i{p,v)T{v^)...T{v2n)) ■ 



(14) 



We are interested in the scaling properties of this quan- 
tity in the regime in which all the separations between 
all the coordinates r, r3, . . . r2„ are of the order of R. For 
p R we can write 

^2,u(r + |,r-|,r3...r2„)~ (|)"'52„(i?) , (15) 

where xi is a yet unknown exponent which in general may 
depend also on n. This exponent will be found below in a 
particular model and will be shown to be n-independent. 
For p very small we can use (0) and (111) to write 

J^2n,l{r + - ■|,r3 . . .r2„) = PaiPa2 ■■■Pai (16) 

oo 

><Ef^^P^^^.t.;'(r,r3,...r2„). 

Finally , in the limit p ^ rj we use the fact that J-2n is 
smooth in p up to p ^ 77 to evaluate the differential op- 
erator as divisions by ry: p'^PW^P ~ (p/?/)^^. Accordingly 
we we have in this limit 



lim J^2n,iir + - ■^,r3 . . .r2„) = 

^1,0 y-ya-, ...Oil { \ I 

Pa,Pa2 • •■Pa,— <-2n,i,0 (l"'I"3,-- •r2n) CX p 



(17) 



Next we want to explore the scaling behavior of J-2n,i 
for values of p in the inertial range rj <^ p <^ L. This we 
cannot do in general. We need to specialize now to a par- 
ticular dynamical model. We choose Kraichnan's model 
of passive advection of a scalar field T(r, t) by a random 
velocity field whose statistics are Gaussian, and whose 
correlation functions are scale invariant in space and 5- 
correlated in time |l^,|l^ . The relevance of the results to 
Navier-Stokes turbulence will be discussed later. For a 
scalar diffusivity k the dissipation field is e(r) = KjVTp 
and the quantities (0) and (2|) are replaced by 



^2„(i?)^(inR)-no)r)^[^2(i?)]"(|) ^ a^) 



if,,(i?) = ((e(R)-e-)(6(0)-e-))^e-2 



R 



(19) 



(v[/(p,r)T(r3)...r(r2„)) 



(13) 



In the present case the scaling exponent of 5*2,1 is C2n = 
nC,2 - <5„. 

It was shown in [p^ that the correlation function !F2n 
solves a particularly simple equation when two of its coor- 
dinates (say ri and V2) are much closer to one other than 
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all the rest. Explicitly, for p small the p dependence of 
this function is governed by the equation 

B{p) J^2n{r+ - ■|,r3 . . .r2„) = $2n-2(r,r3 . • .r2„) . 

(20) 

Here $2n-2(r, . . . r2„) is a homogeneous function with 
scaling exponent C,2n ^ (2- In 3-dimensions the operator 
Bip) is given by 0,01 



B{p) ^ H 



d A-G.d (4-C2) .2 

p^dp' dp 2pC2 



(21) 



Here if is a constant. It has been shown [T^JlJl that the 
leading scaling solution for the p dependence of function 
J-2n is an eigenfunction of the operator B{p) with eigen- 
value and thus can be expanded in spherical harmonics: 



^2n{r + p/2, r - p/2, ra . . . r2„) 

00 I 

E At'(r,r3...r„)pAll,™(p), 



(22) 



1=0 ■ 



where Af^^ (r, r3 . . . r„) is a homogeneous function whose 
scaling exponent is C,2n ~ A . To compute the exponents 
Pi for ^ 7^ we need to find a solution of the homogeneous 
part of ( pO| ) . By a direct substitution of ( p2| ) into the LHS 
of (|^ onefinds A(/3,+3-C2) = (4 - C2)K' + l)/2- Note 
that the LHS of this relation originates from the radial 
part of the operator B, whereas the RHS results from 
the angular part that is proportional to . Solving the 
quadratic equation for /?; we find in 3-dimensions p3|,p^ : 
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C2 - 3 + v/(3-C2)2 + 2/(? + l)(4-C2)l . (23) 



The multipole decomposition of ( p^ , similarly to (^ and 
(^, leads to 



J^2n,i (r + p/2, r - p/2, ra . . . r2„) 

m— — l 



(24) 



In the situations in which all the separations between the 
coordinates r, r3 . . . r„ are of the same order of magnitude 
_R, and i? ^ p ^ r/ we can write 



^2„,;(r + p/2, r - p/2, r3 . . . r2„) cx 



(25) 



Comparing with Eq.(^5|) we identify the exponent xi as 
f3i and write the final form: 



.F2„,,(r+|,r-|,r3...r2„)~ (^-^j ^2«(i?) ■ (26) 



At this point we want to match the solution (26) which 
is valid for p ^ rj with the solution (^ which is valid 



for p <^ rj. This can be done if the solution is varying 
smoothly across 77 without any non-monotonic behavior. 
The rigorous proof of this property is beyond the scope 
of this Letter. It can be demonstrated numerically by 
solving the ordinary differential equation rt2Q). Here we 
simply assume this. Equating ( p5| ) and ( |17| ) for p = rj 
(up to an unknown _R-independent coefficient) we find 



C2n,;,o(r, I"3i ■ • • r2n) 



1 



V \R 



S2n[R) , (27) 



where we remind the reader that R stands for the order of 
magnitude of all the separations between the coordinates 
of C2„,;,o- Comparing with Eq.(|l8|) we conclude that the 
correlation function of L/^o with any even number of T 
fields separated by distances of the order of R depends 
simultaneously on two renromalization scales, rj and L, 
and on the two sets of anomalous exponents Pi and Cn- 

Until now we considered correlations of one new local 
field with a number of T fields. We can also examine a 
cross correlation of two (generally different) local fields. 
Repeating the analysis one finds 



(LLo(r + R)L,,,o(r)) -^77 



-i+V fH 
R 



S^iR) . (28) 



We see that in general such correlations depend on the 
two renormalization scales and on two sets of exponents. 
It is therefore interesting to ask why this phenomenon 
is absent in K^e which is closely related to such cor- 
relation functions. We note that in terms of our lo- 
cal fields the correlation ([l9| ) is given by K^^{R) — 
(Lo,i(r -I- R)Lo,i(r)} as can be checked by substitut- 
ing the definition of the local fields. This is a very special 
case among the correlations of the local fields. The reason 
is that iov I — Q the differential equation (|20|) is inhomo- 
geneous, and it is easy to see by power counting that /5o 
is replaced by C,2- From Eq.([l0|) it follows that in this 
case 



e-M4) . (29) 



RJ 



In the last step we used the fact that by definition e 
— Klim 



ri-raH?) VlV2-?^(ri,r2). SiuCC T{vi,Y2) ^ \yi - 



V2['^ we get e oc nrj' 



C2 



This leads directly to the final 



step in (|29|), in which the renormalization scale 77 from 
the correlator of kLq 1 ^ kV^T^ . The deep reason for 
this is that this is the rate of dissipation of the integral 
of motion in the passive scalar problem. In this sense 
fCee is unusual, and all the generic correlations (^) are 
simultaneously dependent on two renormalization scales. 

One can generate more local fields that will have scal- 
ing properties which may depend on new exponents. In- 
stead of starting with the fusion of two points we can 
fuse three, four or more points. Instead of (|^) we can 
introduce 
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«'3(pi,p2,r) = T(r + pi)T(r + p2)T(r-pi -P2), (30) 

^'4 oc etc. Expanding these fields in Taylor series 
with respect to pi,P2 etc, we can generate new sets of 
local fields that contain derivatives of three, four etc.T- 
fields. Their correlation functions will depend on the 
ultra-violet exponents which appear due to three-point, 
four- point etc. coalescing clusters, and on the infrared 
scaling exponents of six, eight and more-point correlation 
functions. Of course, the actual values of the exponents 
depend on the dynamical model, but the structure of the 
theory is general. To stress this generality we make now 
a few comments about the Navier-Stokes problem. In 
dealing with Navier-Stokes turbulence we need to worry 
from the beginning about Galilean invariance in addition 
to the SO (3) symmetry group. To this aim we will con- 
sider local fields that originate from the fusion of gradient 
fields. The simplest object is 



*f"'(p,r) 



aua(r + p/2) du^{r- p/2) 



dps 



(31) 



From this point on we can proceed following the route 
sketched above for the scalar case. Representing this 
field as a multipole decomposition with respect to the 
direction of p, and considering the Taylor expansion in 
p, we can generate infinitely many local fields. These 
fields have two u fields and as many gradients as we 
want to consider, starting from two. It is interesting to 
notice that in the present case we have two different vec- 
tors, i.e. V and u from which we can form untisymmet- 
ric combinations, like the vorticity = e^fs-fdup/dr^ 
(where Cq^^ is the fully antisymmetric tensor). Conse- 
quently we will have odd as well as even I components 
in this scheme. In addition we have symmetric com- 
binations of velocity derivatives like the strain tensor 
SaB = [dua/drp + duf3/dra]/2. In general the tensor 
( pl| ) has 36 independent components, serving as a ba- 
sis for a 36-dimensional reducible representation of the 
0(3) group (SO(3)-|-inversion). This basis may be de- 
composed into a set of irreducible bases of lower dimen- 
sions. There are two scalar fields, LOaUJa and = Sajjspa 
each of which is a basis for one-dimensional irreducible 
representation with I — 0. The pseudo-vector Sapuj/s is a 
three-dimensional basis for an irreducible representation 
with 1 = 1. There exist three traceless tensor fields each 
of which is a five-dimensional basis belonging to I = 2 
and taking care of 3x5 = 15 components. An example is 



(32) 



In addition we have one 3-rank pseudo tensor correspond- 
ing to / = 3 and one 4-rank tensor corresponding to ^ = 4. 
The last two fields exhaust the remaining 7-1-9 compo- 
nents. As in the scalar case there are fields with all values 
of / which are obtained when more gradients act on our 
field (|l|). Finally we can also start with a higher number 



of fusing gradient fields dua/dpp to generate new sets of 
local fields having three, four and more velocity fields. 
The exploration and utilization of this rich structure is 
beyond the scope of this Letter. It will suffice here to 
state that these local fields will have correlation function 
with anomalous scaling properties that generally depend 
on two renormalization scales and on two sets of anoma- 
lous scaling exponents. The correlator (|^) will again be 
special and ry-independent since it involves the rate of 
dissipation of the integral of motion (energy). Correla- 
tions of fields with I ^ will be generic. For example the 
correlation of VO2 with e = 2vs^ is 



of (r-fR)s^(r) 



(I)' 



(33) 



Our guess is that x is numerically close to p and that y is 
numerically close to 2/3, with an accuracy which is of the 
order of the difference between C2 and its K41 estimate of 
2/3. We stress however that the main point of this Letter 
is not the numerical value of this or that exponent, but 
that normal scaling, which is based on dimensional anal- 
ysis (like K41 for Navier-Stokes turbulence), fails dou- 
bly due to the explicit appearance of two physically im- 
portant scales, the inner and the outer renormalization 
scales. 
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